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ABSTRACT 

Where particle acceleration and plasma heating take place in relation to magnetic reconnection is a 
fundamental question for solar flares. We report analysis of an M7.7 flare on 2012 July 19 observed by 
SDO/ AlA and RHESSI. Bi-directional outflows in forms of plasmoid ejections and contracting cusp- 
shaped loops originate between an erupting flux rope and underlying flare loops at speeds of typically 
200-300 kms~^ up to 1050 kms~^. These outflows are associated with spatially separated double 
coronal X-ray sources with centroid separation decreasing with energy. The highest temperature 
is located near the nonthermal X-ray loop-top source well below the original heights of contracting 
cusps near the inferred reconnection site. These observations suggest that the primary loci of particle 
acceleration and plasma heating are in the reconnection outflow regions, rather than the reconnection 
site itself. In addition, there is an initial ascent of the X-ray and EUV loop-top source prior to its 
recently recognized descent, which we ascribe to the interplay among multiple processes including 
the upward development of reconnection and the downward contractions of reconnected loops. The 
impulsive phase onset is delayed by 10 minutes from the start of the descent, but coincides with the 
rapid speed increases of the upward plasmoids, the individual loop shrinkages, and the overall loop- top 
descent, suggestive of an intimate relation of the energy release rate and reconnection outflow speed. 
Subject headings: acceleration of particles — Sun: flares — Sun: UV radiation — Sun: X-rays, gamma 
rays 



1. INTRODUCTION 

Magnetic reconnection is believed to be the primary 
energy release mechanism during solar flares, but where 
and how the released energy is transformed to heat 
the plasma a nd accelerate particl e s remain unclear (for 
reviews, see iHolman et all 120111: iFletcher et all 120111 : 
lPetrosianll201 2: Ravm ond et al.|[2QT2l ). Evidence of mag- 
netic reconnection and current sheets on the Sun has 
been observed in various situations and wavelengths. 

A major advance in the last decade was the discov- 
ery of a second coronal source above a commonly ob- 
served loop-top source in X-ray s and r adio wavelengths 
jLSui & Holman ^2003"; Sui et al. ^2004"; "Pick et al." ^20051; 
Veronig et al. 2006; Li & Gan 2007; Liu et al. 200i 
2009c; C hen Petrosianll2012l : ISu et al.ll2012HBain et all 
2012: Glesener et al.|[2012f ). Such double sources often 



served in soft X -rays (SXRs; iForbes fc ActonI 119961: 



.eeves et al.ll2008r) , ext reme ultraviolet fEUV: | Li fc Gan 
2006f) . and microwaves (jLi fc Ganll2005l : iReznikova et al 
20101 ) . They were interpreted as contractions of newly re- 



exhibit higher-energy emission being closer to each other, 
indicating higher temperatures or harder spectra of ac- 
celerated electrons in the inner region nearer to the pre- 
sumable magnetic reconnection site. 

Another surprise has been the descent of the loop-top 
X-ray source at typically 10-40 km s~^ early in the im- 
pulsive phase before its common ascent throu gh the de- 
cav p hase (ISui fc Holmanl l2003l: iSui et al.l[200llLiu et all 
[200l[Sh^et al.ll2008[ ). The upper coronal source, how- 
ever, usually keeps ascending all th e time. 

Shrinkages of entire flare loops (jSvestka et al.l 119871 ) 
at speeds on the order of lOkms"^ have been ob- 

^ Lockheed Martin Solar and Astrophysics Laboratory, Build- 
ing 252, 3251 Hanover Street, Palo Alto, CA 94304 

^ W. W. Hansen Experimental Physics Laboratory, Stanford 
University, Stanford, CA 94305 

^ Department of Physics, Stanford University, Stanford, CA 
94305 



connected loops due to magnetic tension as they evolve 
from initially cusp shapes toward more relaxed round 
shapes. 

Often after the impulsive phase, bright loops and dark 
voids seen in SXR or EUV descend onto a flare ar- 
cade from above at greater speeds of typicall y 150 kms~ ^ 
(McKenzie fc Hudson 1999; Savage fc McK enzid [20TTh . 
At even greater heights of a few i?©, similar descend- 
ing loops are seen in white-light coronagra phs, usually 
hours after a coronal ma ss ejection (CME; I Wang et al.l 
119991 : iSheelev et al.l[200l ). These features are also inter- 
preted as contracting post-reconnection loops, but their 
observed speeds are only a small fraction of the expected 
coronal Alfven speed on the order of ~1000kms~^. 

Imaging and Doppler observations have also re- 
vealed bi-directional niagnetic reconnection inflows 
(lYokovama et al.l 120011: iMim^an et al.l [20T0I: iLiu et ail 
2OIOD. out fl ows Cinnesetal. '19971: iKo et all 12003 



Wang et all 120071: Nishizuka et alj 120101: iHara et all 
20111: IWatanabe et al.l I2012D . or both (iLin et al.l l2' 



Takasao et all 120121 : iSavage et ani2012f ). 



The physics behind descending X-ray loop-top sources 
and their relationship with slow loop shrinkages, 
fast supra-arcade loop contractions, and reconnec- 
tion outflows remain unclear, although some mod- 
els have been proposed (e.g., iSomov fc Kosugil 119971 : 
iKarlickv fc Kosugil l2004f ). To fill this gap, we present 
here observations of a recent eruptive M7.7 flare from 
the Reuven Ra maty High Energ y Solar Spectroscopic Im- 
ager ( RHESSI: \Lin et al.lIMM and Solar Dynamics Ob- 
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servatory Atmosp heric Imaging Assembly {SDO / AlA; 
iLemen et alJl2Ql"2l ). We find in this flare all the above 
interrelated phenomena, which can be understood in a 
unifled picture as contractions of post-reconnection loops 
modulated by the interplay between energy release and 
cooling. The observed highest speed of ~1000kms~^ of 
loop contractions is comparable to the expected Alfven 
speed. The maximum temperature and nonthermal loop- 
top emission being away from the inferred reconnection 
site suggest that primary heating and particle accelera- 
tion take place in the outflow regions, rather than the 
reconnection site itself. 

After an observational overview in Section [21 we 
present motions of the overall X-ray and EUV emission 
in Section [3l We examine bi-directional outflows in forms 
of plasmoids and contracting loops in Section HJ In Sec- 
tion O we analyze the spatial distribution of energy and 
temperature dependent emission, including double coro- 
nal X-ray sources. We conclude in Section [6l followed 
by two appendixes on supplementary AIA and STEREO 
observations. 



TABLE 1 
Event Time Line (2012 July 18-19) 

22:18, 07/18 Peak of the earlier C4.5 flare 

04:17, 07/19 Onset of the M7.7 flare and initial ascent 

of the overall X-ray and EUV loop-top source 
05:02 05:07 Onset of overall X-ray and EUV loop-top descent 

and transition from slow to fast rise of the flux rope CME 
05:15 Max. velocity (1050 km s"""^) of upward ejections 

05:16 Max. velocity ( — 58kms~"'^) of downward loop shrinkages 

05:15 05:20 Max. velocity (-7 to -23kms~^) of overall 

X-ray and EUV loop-top descent 
05:16-05:43 Flare impulsive phase, hard X-ray burst 
05:21-05:31 Min. height of overall X-ray and EUV loop-top descent 

and onset of the second ascent 
04:17 05:20 Upward ejections; downward, low-altitude fast contractions 
05:20 16:00 Downward, high-altitude fast loop contractions 
04:17 16:00 Downward slow loop shrinkages 

06:52 Max. velocity (— 918kms~^) of downward fast contractions 



2. OVERVIEW OF OBSERVATIONS 

The event under study was an M7.7 flare that occurred 
at -04:17 UT on 2012 July 19 in NOAA active region 
(AR) 11520 on the southwest limb. It was well observed 
by RHESSI and SDO/AIA, but it was not detected by 
Fermi and its impulsive phase was missed by the X-ray 
Telescope on Hinode. Table [T] summarizes the event time 
hue that will be discussed in detail. 

Figure [T] shows the history of the flare emission. The 
GOES 1-8 A flux peaks at 05:58 UT followed by a slow 
decay lasting almost one day. We define the interval of 
05:16-05:43 UT as the impulsive phase, as marked by the 
two vertical dashed lines, which starts at the sudden rise 
of the RHESSI 25-50 keV flux and ends (during RHESSI 
night) when the time derivative of the GOES 1-8 A flux 
dro ps to its level at the impulsive phase onset, assuming 
the iNeupertI (| 19681 ) effect at work. We call the intervals 
before and after the impulsive phase the pre-impulsive 
and decay phases. RHESSI hdiS good coverage except for 
the late impulsive and early decay phases. 

Figure [2] and its associated Movie A show AIA images 
of the event. An earlier C4.5 flare occurred in the same 
location, peaking at 22:18 UT on 2012 July 18 (see Fig- 
ure [TJa)). This is a confined flare that produces a hot 
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Fig. 1. — History of X-ray and EUV flare emission, (a) GOES 
SXR fluxes showing the preceding C4.5 flare and the main M7.7 
flare, (b) SDO/ AIA intensity of the loop-top region (at projected 
height /iref = 60^^ on Cut 0, as shown in Figure I13|) normahzed 
by its initial value (base ratio). Color-coded for AIA channels, the 
curves are vertically shifted to avoid overlap and arranged from 
the top to bottom approximately in the temporal order of their 
initial response to the flare (in progressively cooler channels), (c) 
RHESSI count rates in colored solid lines and GOES 1-8 A flux 
and its time derivative in black dotted lines, arbitrarily shifted ver- 
tically. The two vertical dashed lines indicate the impulsive phase, 
(d) Temperature and (e) emission measure of the flare plasma in- 
ferred from RHESSI (orange/red) and GOES (black) spectral flts, 
together with the power-law index 7 of the nonthermal component 
shown in blue in (d). 

flux rope failing to erupt and cusp-shaped flare loops 
underneath it (Figure [2fa)). This configuration then 
gradually evolves for hours and finally becomes unsta- 
ble, initiating the later, eruptive M7.7 flare, when the 
flux rope is expelled as a fast CME of >1000kms~^. 
The flux rope in this two - stage eruption was reported by 
iPatsourakos et all (j2Q13r ). 

As shown in Figure [2fi), the trailing edge of the CME 
displays a clear "V-shape", which, together with the 
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underlying "inverted- V shape" of cusp-like flare loops, 
suggests two Y-type null points with a vertical current 
sheet formed in between, as predicted in t he clas sical pic- 
ture of eruptive flares (ICar michael 1963: iStiirro ck 1966; 
lHiravamalll974l : iKopp fc Pneumanlll976[ ). Not predicted 
in that picture is the initial upward growth of the cusp 
followed by its rapid downward shrinkage around the 
early impulsive phase, prior to another, expected up- 
ward growth through the decay phase in a commonly 
observed candle-flame shape (Figures [2]^a)-(l)). The ini- 
tial growth and shrinkage are accompanied by the grad- 
ual rise and impulsive eruption of the overlying flux 
rope, respectively. Equally interesting are high-speed bi- 
directional outflows involving upward-moving plasmoids 
and downward-contracting pointed cusps (Figure [21 bot- 
tom; Figure [5]). We examine these and related features 
observed by RHES S I dind AIA in next several sections. 

3. OVERALL X-RAY AND EUV SOURCE MOTIONS 

We flrst follow the evolution of the morphologies and 
positions of overall X-ray and EUV sources. 

3.1. X-ray Source Morphology 

We reconstructed RHESSI X-ray images in energy 
bands from 3 to 50 keV using the CLEAN algorithm and 
detectors 3-9. Depending on the count rate, we chose 
variable integration time ranging from 20 s during the 
impulsive phase to 4 minutes during the decay phase. 

Figure [3] and its accompanying Movie show examples 
of RHESSI contours overlaid on AIA images. There is 
a persistent loop-top source at 6-10 keV (green) near 
the apex of the cusp-shaped EUV loops. Accompanying 
the evolution of the EUV loops mentioned above, this 
loop-top source undergoes a gradual ascent followed by a 
descent and then another ascent. This can be best seen 
in the last panel showing the temporal migration of the 
emission centroid obtained from a contour at 70% of the 
maximum of each image. Early descen ts of loop-top X - 
ray sources have been recognized (e.g., iSui et al] 120041 ). 
but this is the flrst time that an evident preceding ascent 
is observed. 

Early in the event there is an additional, weaker coro- 
nal X-ray source (panel (b)) located in the lower por- 
tion of the overlying flux rope. It later falls below 
detection with the eruption of the flux rope. Dur- 
ing the impulsive phase (panel (e)) at higher energies 
(25-50 keV), th ere are double footpoint sources and 
a Masuda-type (|Masuda et all 119941 : iNitta et al.l I2010L 
e.g.,) coronal source located 20^' = 15 Mm above the SXR 
(6-10 keV) loop-top. This separation is twice that of the 
Masuda case. 

The earlier C4.5 flare (panel (a)), though more com- 
pact, exhibits surprisingly similar emission, including the 
additional X-ray source within the flux rope and the 
Masuda-type hard X-ray (HXR) source. This suggests 
that the end state of the flrst, conflned flare of a failed 
eruption serves as the initial state of the second, erup- 
tive flare, making them homologous flares. Note that 
the weaker, southern footpoint emission is (partially) oc- 
culted by the Hmb, especially for the flrst flare. 

3.2. X-ray and EUV Height-time History 

To track various moving features, we placed eight cuts 
of 10^' wide centered on the limb, as shown in Figure [2jf). 



Cut is positioned along the presumable current sheet 
between the initial flare cusp and erupting flux rope, and 
is used as the flducial line for obtaining projected heights 
up to ~08:00 UT. Subsequent cuts are evenly spaced by 
3° to capture the cusp at different stages as it gradual 
turns toward the south. For each cut, we averaged pixels 
of an AIA image sequence within it in the perpendicular 
direction to obtain a space-time plot. 

Figure m^c) shows, for example, a base-ratio (i.e., nor- 
malized by the initial intensity proflle) space-time plot 
from Cut at 131 A. It shows the eruption of the flux 
rope and the three-stage development (upward, down- 
ward, and upward again) of the underlying flare cusp. 
The peak emission at each time, marked by the small 
purple symbols, evidently exhibits this motion. Space- 
time plots of other AIA channels are shown in Figure [13] 
and described in Appendix A. 

Figure Hl^e) shows the height-time history of the loop- 
top emission centroids in four RHESSI X-ray bands from 
3 to 25 keV and in all AIA EUV channels. (From now 
on, we refer to both X-ray centroids and EUV peaks as 
centroids.) We flnd that all loop- top centroids follow 
the same three- stage motions. The AIA 131 and 94 A 
channels best exhibit this trend continuously, while in 
cooler channels, especially 211 and 171 A, the initial as- 
cent and subsequent descent become obscure. Note that 
the 3-6 keV centroid (red diamonds) cannot be identi- 
fled since the impulsive phase because the RHESSI thin 
attenuator moves in and raises the energy threshold to 6 
keV. 

We also notice a clear energy dispersion that the loop- 
top emissions at higher temperatures or photon energies 
are located at greater heights, with all X-ray centroids 
situated above EUV centroids. For example, at its great- 
est height prior to the impulsive phase, the 16-25 keV 
centroid is 35'' = 26 Mm above its 335 A counterpart. 
This energy dis persion is consisten t with previous obser- 
vations (e.g., iVeronig et al.l [2QQ6I ) and in line with the 
expected picture that higher loops are newly energized 
and thus are hotter and/or host nonthermal electrons of 
harder spectra, while lower loops are previously energized 
and have undergone cooling. We suggest that the so- 
called above-the-loop-top Masuda-type sources, includ- 
ing the 25-50 keV coronal source in Figure [3fd), could 
be extreme cases of this general trend. 

In Table [21 we list the initial heights of X-ray and 
EUV loop-top centroids around 04:35 UT, together with 
the maximum and minimum heights during their de- 
scents and the percentage height reductions. In gen- 
eral, the above noted energy dispersion persists and the 
10-20% descents ar e in agreement with earlier reports 
(|Veronig et al.ll2006[ ). An exception is the 193 and 304 A 
channels because of their response to both hot and cool 
emissions (see Appendix A). There is a weak trend that 
higher energy X-ray descents start (at the maximum 
heights) later, but still within 5 minutes during 05:02- 
05:07 UT, and X-ray descents end (at the minimum 
heights) earlier than EUV descents by up to 10 minutes. 

We measured the average velocities of the loop-top cen- 
troids by fltting a piecewise linear function to the height- 
time data during their initial ascent, subsequent descent, 
and second ascent stages. As listed on the right hand side 
of Table [21 the result shows a general trend of higher ve- 
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Fig. 2. — AIA images of the event, rotated to the solar west up. Top: 94 A image sequence of the earher C4.5 flare and the M7.7 flare 
under study (see the associated onhne Movie A for composite 94 (red) and 335 A (green) images). The brackets in (c) and (e) mark the 
enlarged flelds of view (FOV) of the middle and bottom panels, respectively. The numbered dotted lines in (f) mark the cuts for obtaining 
space-time plots presented in this paper. Middle: Detailed 131 A image sequence showing the upward growth of the flare cusp, followed 
by its rapid shrinkage around the early impulsive phase. The horizontal dotted line marks the highest position of the early cusp. Bottom: 
131 A image sequence showing an example of pointed cusps contracting from a ray-like structure, presumably a reconnecting current sheet 
(see the associated Movie B). The flrst panel is an original image and the rest are running difference images. The slanted dashed line 
indicates an average contraction velocity of — ISlkms"-*^. 



locities at higher energies. The X-ray descents of —11 
to — 15kms~^ are about twice faster than the EUV de- 
scents except for 193 A because of its dual temperature 
response. There is also a gradual decrease in velocity 
with time during the second ascent stage through the 
decay phase. These trends generally a^ ree with previous 
observations (|Liu et al. 2QQ1 : IVeronig et al. 2QQ6i) . 

To follow the temporal variations of the loop-top ve- 
locities, we took time derivatives of spline fits to the cen- 
troid heights of selected channels. The resulting veloc- 
ities vs. time are shown in Figure HJ^c) for the 6-10 & 
10-16 keV and 131 & 94 A channels. The outstand- 
ing feature is that the maximum descent velocities (also 
tabulated in Table [2]) in the range of —7 to — 23kms~^ 
are all attained during 05:15-05:20 UT, near the onset 
of the impulsive phase at 05:16 UT marked by a ver- 



tical dashed line. This time is delayed from the initial 
descent by about 10 minutes. This indicates that the 
loop-top descent velocity, rather than its height, is more 
intimately correlated with the energy release rate. 

We also show in Figure lU^e) the centroid heights of 
the upper coronal X-ray source located in the lower por- 
tion of the flux rope. A parabolic fit indicates a fi- 
nal velocity of 19kms~^ at 04:50 UT, comparab l e to or 
slower than those in other events (|Liu et al.ll2008L l2009cl : 
iSui fc Holman|[2"003f ) . It is also 50% slower than the mid- 
dle portion of the then accelerating flux rope, as indi- 
cated by the parabolic fit in Figure [THTe). which reaches 
153 ± 3kms-^ later at 05:07 UT near the edge of AIA's 
FOV. 

The conjugate footpoints at 25-50 keV during the im- 
pulsive phase move away from each other at an average 



Particle Acceleration in Reconnection Outflows 



5 




220 230 240 250 260 
-Solor Y (orcsec) 



(b) 04:42:35 (c) 05:06:35 (d) 05:15:11 (e) 05:25:01 



(f) 06:15:49 (g) X- ray LT Centroids 




200 220 240 260 
-Solor Y (orcsec) 



200 220 240 260 
-Solor Y (orcsec) 



200 220 240 260 
-Solor Y (orcsec) 



200 220 240 260 
-Solor Y (orcsec) 



200 220 240 260 
-Solor Y (orcsec) 



225 230 235 240 



Fig. 3. — (a)-(f) Evolution of RHESSI 6-10 keV sources (green contours at 31% and 90% of the maximum) overlaid on concurrent AIA 
131 A images (see the accompanying Movie). The blue contours show 25-50 keV emission, for both the C4.5 and M7.7 flares, at the 
footpoint(s) and the Masuda-type coronal source 2Q" above the SXR and EUV loops, (g) Color-coded temporal migration of the 6-10 keV 
loop-top centroid as marked by the plus sign on the left. The red arrows indicate an initial ascent followed by a descent and then another 
ascent. The dotted line marks the fiducial line (Cut 0) for obtaining projected heights shown in Figure [H The first and last panels are 
enlarged from the bracketed regions shown in their adjacent panels. 



TABLE 2 

Projected Heights and Velocities of RHESSI and AIA Loop-top (LT) Centroids 

Channels Height (arcsec^O. 737 Mm)/Time (UT) Velocity ( km s~^)/Time (UT) ~ 

Initial Max Mm Descent % Initial Ascent Descent 2nd Ascent 
(^04:35) Mean Max 6:00 6:00 7:00 7:00 8:00 
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Note. — Centroid height uncertainties are <2" for RHESSI and <l" 

velocity of 13kms~^ (see Figure Hl^e)), almost identi- 
cal to the velocity of the simultaneous loop-top ascent 
(see T able [2]). T his is consistent with previous observa- 
tions (|Liu et alJ[2QQ4) and the classical picture of erup- 
tive flares during the arcade growth phase. 

4. BLDIRECTIONAL PLASMA OUTFLOWS 

We now turn our attention to the bi-directional plasma 
outflows observed by AIA. As shown in Figure [5] and its 
accompanying Movie, emission features move both up- 
ward and downward from above the cusp-shaped flare 
loops. The upward outflows, observed only when the flux 
rope rises toward its eruption, are mainly bright blobs 
(interpreted as plasmoids, top row), while the down- 
ward outflows are primarily in the form of retracting 
cusp-shaped loops (bottom row) observed throughout 
the event well into the decay phase. Among the down- 
ward retracting loops, those seen at low altitudes tend to 
have lower speeds (<60kms~^) and gradually decelerate 
when approaching and piling up onto the apex of the flare 
arcade, while those originating from high altitudes tend 
to have higher speeds (>100kms~^) and decelerate more 
rapidly. According to these observational distinctions, 
we call the former slow loop shrinkages and the latter 
fast loop contractions^ although they might share a com- 
mon physical origin as discussed in Section [6l The former 
are likely analogous to shrinkages observed at SXR and 
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for AIA. Velocity uncertainties are on the order of 10% 15%. 

other wavelengths (|Svestka et al.l 119871 ) , while the lat- 
ter, especially those occurring after the flare arcade has 
formed, are likely so-called supra-arcade downflow loops 
(|Savage k McKenzidlMTl ). 

To track these moving features, we obtained space-time 
plots from the cuts defined in Figure [2ff) using the 131 A 
channel because it provides the best coverage of these fea- 
tures. We then applied running ratio, namely, dividing 
the intensity profile at each time by its previous neigh- 
bor, which highlights moving features as intensity tracks 
(see, e.g.. Figure Hl^d)). We exhaustively identified such 
tracks of more than 10'^ long corresponding to unique 
moving features. If one feature is captured by multiple 
cut s, we included only t he most complete track. Follow- 
ing l^^rre^^F^l] (j2Qlll ). we fitted the projected height 
h of each track by a function of time t, 

hit) =ho^ VTt + aor\e-'^^ - 1), (1) 

which gives the instantaneous velocity 

v{t) = VT - aore-'^\ (2) 

and acceleration 

a{t) = aoe-^/^, (3) 

where ho and ao are the initial height and acceleration, 
vt is the terminal velocity, and r is the e-folding decay 
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Fig. 4. — X-ray and EUV source motions, (a) and (b) Sub- 
sets of light curves from Figures [TJb) and (c). (c) and (d) Base 
and running ratio space-time plots of the AIA 131 A channel from 
Cut 0, overlaid with selected centroid positions from (e). (e) Pro- 
jected heights of the emission centroids of the upper and lower 
coronal sources from RHESSI (large symbols, greater heights) and 
of the loop-top peaks from AIA (small symbols, lower heights). 
Also shown is the 25-50 keV footpoint separation (purple squares) 
shifted upward by 30^^ (f) Velocity of centroids from spline fits to 
selected channels in (e). The blue symbols are the initial velocities 
of individual but rescaled by a factor of 0.5. 
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Fig. 5. — AIA 131 A images (see the accompanying Movie) show- 
ing upward moving blobs (top, labeled A and B) and V-shapes 
(middle) and downward contracting cusps (bottom). Their space- 
time fits are shown in Figure [Tfa). The green contours show two 
distinct 8-10 keV RHESSI sources from Figure [Tolc). The first row 
and column are original images and the rest are running difference 
images. 

time that characterizes the observed decreasing acceler- 
ation (or deceleration) with time. We present below the 
kinematics of the features categorized above and focus 
on their initial velocities vq = vt — clqt and accelerations 
ao. 

4.1. Slow Downward Loop Shrinkages 

Figure [6] shows kinematic measurements of the slow 
loop shrinkages. These loops traverse a hot (up to 
T ~ 30 MK; see Figure [TT|) loop-top region, which is 
bright only in the hottest channels of 193, 131, and 94 A, 
and descend toward the apex of the flare arcade seen in 
cooler channels, where they fade below detection (see 
Figure [T3|) . They generally start at a low velocity with a 
median of — 17kms~^ (+/~ upward/downward) and 
decelerate to typically — 5kms~^ within 0.5-2.5 h. Their 
median initial deceleration is 0.014 km s~^. 

Such shrinkages are persistent throughout the flare, 
not only during the decay phase when the overall 
loop-to p emission develops u pward as previously ob- 
served (jForbes fc Actonlll996l ). but also during the pre- 
impulsive and impulsive phases when the loop-top un- 
dergoes its initial ascent and subsequent descent. 50 
shrinkages are identified during 04:00-16:00 UT with an 
average occurrence rate of once per 7.7 minutes up to 
10:00 UT but more frequent around the impulsive phase 
(Figure [6i;b)). 

Regardless of the direction of motion of the loop-top 
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Fig. 6. — Kinematic measurements of slow loop shrinkages, (a) Running ratio space-time plot enlarged for the boxed region in (b), 
overlaid with examples of fits to shrinkage tracks labeled with their initial velocities, (b) Base ratio 131 A space-time plot from Cut 
overlaid with all 50 fits performed, (c) Initial velocities vq of the fitted tracks in (b) as a function of time, (d) Histograms of vq (solid) and 
the final velocity Vf (dotted) measured at the end of each track, with their statistical properties (median, mode, mean, standard deviation, 
and total number of fits) listed in (c). (e) and (f) same as (c) and (d), but for the initial acceleration ao- The black dotted lines in (c) and 
(e) are the time derivative (inverted) and original GOES 1-8 A fiux. The two vertical dashed lines mark the impulsive phase. 



emission, the initial height of loop shrinkage generally in- 
creases with time. The initial velocity vq exhibits more 
variability, as shown in Figure Efc), and appears to be 
positively correlated with the flare energy release rate in- 
dicated by the time derivative of the GOES flux (black 
dotted line). In particular, near the onset of the impul- 
sive phase at 05:16 UT (vertical dashed line), '^o increases 
rapidly and reaches a maximum of — 58kms~^. There is 
a similar temporal correlation with the initial deceler- 
ation ao (panel (d)) that also generally correlates with 

Vq. 

4.2. Fast Downward Loop Contractions 

Figure [7] shows, in the same form as Figure [6l kine- 
matic measurements of fast downward loop contractions 
as well as upward plasmoid ejections that will be ex- 
amined in Section 14.31 We color-code fits to tracks of 
different categories: upward ejections in red and down- 
ward contractions in cooler colors. For the latter, we use 
green for early contractions that occur simultaneously 
with upward ejections and originate from lower heights 
h < 100'', and blue for later contractions that occur with- 
out observed upward counterparts (possibly out of AIA's 
FOV, after the flux rope has erupted) and originate from 
greater heights h > 100''. We identified 29 green tracks 
and 174 blue tracks during 04:00-16:00 UT. 



As shown in Figure [^c), the fast downward contrac- 
tions, especially those later blue-colored ones, start well 
above the hot loop-top region and travel into it with de- 
celeration for some distance before fading below detec- 
tion. Their final heights are generally above the original 
heights of slow loop shrinkages (see Figure [6l^b)). 

In original AIA images, the fast contracting fea- 
tures are usually bright, cusp-shaped loops (see Fig- 
ure [21 bottom), but occas i onally dark, tadp o le-like 
voids (McKenzie fc HudsonI 119991 : llnnes et"a U I200I 
Cassak & Drake 2011). In running ratio space-time plots 
(Figure [71 top), the former each produce a bright track 
followed by a recovering dark track, while the latter pro- 
duce tracks of reversed order. We treat them all as con- 
tracting loops and do not distinguish their differences 
that are beyond our scope. 

As shown in Figure [Tj^c), the early, green-colored con- 
tractions have short ranges of 10"-20", because of the 
compactness of the flare at this time, while the later, 
blue-colored contractions travel great distances up to 
200" = 150 Mm. However, their statistical medians, 
modes, and means of the initial velocities are quite sim- 
ilar within 200-300 kms-^ (Figures [Tj^d) and (e); Ta- 
ble [3]) , suggestive of a common origin of these contrac- 
tions possibly as reconnection outflows. The initial decel- 
erations of the early contractions are, on the other hand. 
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Fig. 7. — Same as Figure [6] but for fast upward plasmoid ejections and downward loop contractions, (a) and (b) Running ratio space- 
time plots obtained from Cuts and 2, enlarged for the two boxed regions in (c). Overlaid with examples of color-coded fits, (a) shows 
simultaneous upward ejections (red dotted) and downward contractions of early, low- lying loops (green dotted), and (b) shows downward 
contractions of later, high-lying loops (blue dotted) from well above the flare arcade. Example images corresponding to the flts in (a) 
and (b) are shown in Figures [5] and [2] (bottom) , respectively, (c) Base ratio 131 A space-time plot from Cut 2 overlaid with all flts to 
these identifled features. The thick green dashed line is a spline fit to the initial heights of bi-directional outfiows, as an estimate to the 
position of reconnection. (d)-(g) Initial velocities vq and accelerations ao of the color-coded fits in (c) as a function of time, together 
with their histograms on the right. The blue dotted line in (e) is for the final velocities of the late, high-lying loop contractions. The 
short, vertical stripes in (a) w ith a periodicity of ~1 minute between 04:45-05:15 UT are possibly due to quasi-periodic fiare pulsations 
(jNakariakov Melnikovll2009l ), not to be confused with the tracks of steep slopes of plasma outfiows. 
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Kinematic Statistics of Bi-directional Outflows 
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0.848 
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Note. — Abbreviations: Med. for Median, Mod. for Mode. ('^) Average occurrence interval up to 10:00 UT. 
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about twice higher than those of the later ones, with a 
median of 2.022 vs. 0.848 kms-^ (Figures [T^f) and (g)). 

The later, high- altitude contractions have a mode 
velocity — 158kms~^ that is similar to the me- 
dian velocity — 15 0kms~^ of supra-arcade do wnflows 
(|Savage fc McKenzie 20 ij; IWarren et all l2011f ) . How- 
ever, their median and mean velocities, —242 and 
— 282kms~^, are nearly twice greater. This is due to 
the high velocity tail with 36 out of 174 or 21% con- 
tractions being faster than — 400kms~^ (hatched area 
in Figure [Tl^e)), reaching a maximum of — 918kms~^. 
We ascribe this difference to AIA's improved capabilities 
that allow us to detect such contractions in their early 
stages when their velocities are high and emissions are 
weak. 

The initial velocities and decelerations of the earlier 
and later fast contractions, when taken together, gener- 
ally increase with time up to ~06:10 UT and then de- 
crease. This maximum time is delayed by almost 1 h 
from that of the slow loop shrinkages and the onset of 
the HXR burst represented by the GOES flux deriva- 
tive shown in Figure [6{c). The temporal variations here 
seem to be more closely correlated with the GOES flux 
itself. The late-phase decreases of the initial velocities 
and decelerations are physically reasonable as energy re- 
lease subsides, but there is an observational bias for those 
tracks starting at the edge of AIA's FOV which may have 
decelerated before they are first detected. 

Figure [8] shows the distributions of fractional height 
reductions of the downward loop shrinkages and con- 
tractions, defined as {h^ — hf)/ho^ where ho and hf are 
the initial and final projected heights. The slow shrink- 
ages (light blue) and later, high-altitude fast contrac- 
tions (blue) have similar distributions with medians of 
0.3 and 0.26, respec tively, which are comparable to those 
observed in SXRs (jForbes fc Acton] 119961: iReeves et al.l 
I2008D and predicted theoreticallv ( jLinll2004f ). The earlv. 
low-altitude fast contractions (green) have a median of 
twice smaller, but comparable to the 10-20% descent of 
the overall loop-top emission (Table Ej). Meanwhile, the 
later fast contractions have a sizable fraction of 16% (27 
out of 174) with height reductions more than 0.4 (up to 
0.7), which is larger than predicted. Note that these ob- 
served values are lower limits, because loops, especially 
those detected at the edge of AIA's FOV, could con- 
tinue contraction before and after their observed inter- 
vals when their emission remains below detection. 

4.3. Fast Upward Plasmoid Ejections 

As shown in Figure [71 upward plasmoid ejections 
appear as steep, elongated tracks (red fits). Com- 
pared with the downward loop contractions (green fits) 
at the same time, they travel greater distances (20''- 
80'' vs. 10"-20"), with larger initial velocities (median: 
322 vs. — 208kms~^) but smaller decelerations (median: 
—0.018 vs. 2.022 kms~^). (The opposite signs of the ac- 
celeration and velocity here indicate deceleration.) Some 
upward ejections even experience acceleration. Such dif- 
ferences, especially the ratio of their median velocities of 
208/322 = 2/3, are consistent with previous observations 
and si mulations of oppositely directed r econnection out- 
flows (iTakasao et al. 2012; Barta et al. "2008; "S hen et al.l 
120111 : iMurphv et ah .2012: .Karpen et al.. .2012.) . These 
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Fig. 8. — Histograms of fractional height reductions of slow 
shrinkages (light blue) and early, low-altitude (green) and later, 
high- altitude (blue) fast contractions. The hatched area indicates 
16% of the later fast contractions with reductions greater than 0.4. 
The cyan triangle marks the descent fraction of the 10-16 keV 
loop-top from Tabled 

differences are likely due to different environments: the 
upflow runs into a region of (partially) open field lines 
with low plasma density, while the downflow impinges 
on high density, closed flare loops that can cause strong 
deceleration. 

The starting heights of these upward ejections and 
their downward counterparts (Figure [71(c)) both increase 
with time. This suggests an upward development of the 
recon nection site situa ted between the opposite outflows 
(e.g. JShen et al.|[2Qllf ). We estimated the height of the 
reconnection site from a spline flt, shown as the green 
dashed line in Figures [7j a)- (c), to the initial heights of 
the bi-directional fast outflows. Its initial velocity of a 
few kms~^ is similar to that of the initial loop- top as- 
cent, suggesting the upward development of reconnection 
being its underlying mechanism. This estimated recon- 
nection site leaves the AIA FOV at about 06:40 UT. 

During 04:15-05:25 UT, we identifled 25 upward ejec- 
tions (red) and 29 downward, early fast contractions 
(green), at average occurrence rates of once every 2.4 and 
2.2 minutes, respectively. These rates are comparable to 
that of the later fast contractions (blue, up to 10:00 UT) 
of onc e every 2.0 minute s and to those in MHD simula- 
tions (jBarta et al.ll2008l )FI The persistence throughout 
the entire flare of such contrac tions also agrees w ith the 
statistical result from Yohkoh (jKhan et al.ll2007[ ). 

The initial velocity of upward ejections, as shown in 
Figure [7Jd), increases with time and seems to be cor- 
related with the height of the overlying flux rope as it 
evolves into eruption. In particular, there is a rapid ve- 
locity increase at the onset of the impulsive phase, reach- 
ing 1050 km s~^. T his a s pect of this event is indepen- 
dently studied by R. [03 (|2013[ ). 

5. SPATIAL DISTRIBUTION OF ENERGY AND 
TEMPERATURE DEPENDENT EMISSION 

5.1. RHESSI X-ray Spectra 



^ The ~2 minute intervals of fast contractions are compara- 
ble to the typical periodicities of recentl y detected quasi- p eriodic 
fast-mode m agnetoso nic wave trains (L iu et all I2011bl . I2012d : 
lOfman et al.l 2011: She n Liul [2012') that are correlated with flare 
pulsations. This may indicate that those waves are triggered by 
the energy release episodes associated with these contractions. 
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To infer the overall properties of nonthermal parti- 
cles and thermal plasma, we analyzed spatially inte- 
grated RHESSI spectra of indiv idual dete c tors fo llow- 
ing the procedures deta iled in iLiu et all ()2QQ8l ) and 
iMilligan fc DennisI <^m). 

Figure [9] shows example spectra in three phases of 
the flare. (1) During the pre-impulsive phase (04:59:00- 
05:01:32 UT), we fitted the spectrum with two isother- 
mal functions combined: a warm component (red dot- 
ted) of temperature T = 13 MK and emission measure 
EM = 0.37 X 10^^ cm~^ plus a hot co mponent (ffleen dot- 
ted) of 30MK (also called super-hot. ICaspi fc Linll2010[ ) 
and a 40 times smaller EM. If we were to replace the hot 
component with a power-law ex ^ we would have a 
steep slope of 7 = 7.6. (2) During the impulsive phase 
(05:25-05:26 UT), we used an isothermal (T = 22 MK, 
EM = 2.8 X 10^^ cm-^) plus power-law (blue dotted, 
7 = 3.8) model. (3) During the decay phase, the data 
can be fitted with one isothermal function of a slightly 
reduced temperature of 18 MK but twice higher EM of 
5.5 X lO'^^cm"^. Note that at £;>25 keV, the power law 
in the impulsive phase is more than an order of mag- 
nitude higher than the thermal component. This indi- 
cates that the 25-50 keV Masuda-type loop-top source 
shown in Figure [3] is primarily nonther mal and likely a 
particle acceleration region itself (e.g., ILiu et al.l [20081 : 
iKrucker et aIll20Tol ). 

RHESSI X-ray Spectra 
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Fig. 9. — RHESSI spectra at three selected times before, during, 
and after the impulsive phase. The lower two spectra are verti- 
cally shifted by decades to avoid overlap. Colored dotted lines are 
isothermal and power-law fits whose parameters are listed on the 
top. The black solid lines are the total fits. 

We applied such spectral fits throughout the flare and 
the temporal variations of the fitting parameters are 
shown in FiguresITJd) and (e). In general, the warm ther- 
mal component is persistent in time, and both its T and 
EM (orange plus signs) increase through the HXR burst 
followed by a more gradual decrease. A similar trend 
is present for the hot thermal component (red crosses) 
during the pre-impulsive phase. A thermal fit to GOES 



data gives a somewhat lower T but higher EM (black 
dotted lines), because of its well-known preferential re- 
sponse to relatively cooler plasma. The power-law com- 
ponent (blue) displays a common hardening trend before 
the RHESSI night data gap. 

5.2. Energy Dependence of Double Coronal Sources 



In Section [321 we have examined the energy-dependent 
height distribution of the loop-top X-ray source (lower 
coronal source). Here we extend this analysis to the up- 
per coronal source during the pre-impulsive phase and 
compare the two sources together. 

As shown in Figures [TOT a) and (b), the higher-energy 
emission (12-15 keV, blue contours) of the lower source 
is located at higher altitudes than lower-energy emission 
(3-8 keV, green contours), while the upper source has an 
opposite trend. That is, the two sources are closer to each 
other at higher energies, as can be better seen in panel 
(d) from the linear fits to their centroid heights as a func- 
tion of logarithmic energy. This implies higher tempera- 
tures of the thermal plasma and/or harder spectra of the 
nonthermal particles within the inner re gions. Similar 
RHESSI observ ations have been reported (jSui fc HolmanI 
I2QQ3I : iLiu et aL 2008, 2009c) and interpreted as evidence 
of magnetic reconnection and associated energy release 
being located between the double coronal sources. 

At the highest energies (^10 keV), this trend seems to 
be reversed, especially for the u pper sourc e. In previ- 
ous observations (e.g.. Figure 4 in lLiu et"al]l2008f ). more 
obvious reversals occurred at slightly higher energies of 
15-20 keV and were interpreted as a transition from the 
low-energy thermal regime to the high-energy nonther- 
mal regime, in which greater stopping distances of higher 
energy electrons can cause higher energy bremsstrahlung 
X-rays being located farther away. We further suggest 
that this reversal may be an indication of the two X- 
ray sources being spatially separated^ instead of merging 
together at the highest energies; so are the locations of 
primary heating and particle acceleration from magnetic 
reconnection. This is supported by the separate temper- 
ature peaks at the two sources as shown in Figures [TlTa) 
and (h) and discussed in Section [Ql 

5.3. AIA Temperature Maps 

To infer the spatial distributions of temperature T and 
emission measure (EM), we employed a forwar d-fitting 
algorithm using AIA filter ratios (lAs chwand en et al 



201: 



for a different approach, see iBattaglia fc Kontar 



iBatt 



20121 ). This algorithm assumes a Gaussian differential 
emission measure that has a peak EM at temperature 
T. Such a model provides the simplest description of the 
EM weighted temperature distribution along the line of 
sight. 

Figure [11] (left) shows examples of T and EM maps. 
Early in the event at 04:49 UT, the highest temperature 
of logT = 6.95 or T = 8.9 MK is located in the ris- 
ing flux rope and the underlying flare cusp, but EM is 
highly concentrated in the latter. The temperature then 
increases with time especially in the flare cusp, reach- 
ing logT = 7.45 or T = 28.2 MK at 05:09 UT, close to 
the 35 MK of the hot RHESSI component at this time 
(see Figure [l](d), red crosses). In the decay phase, the 
highest temperatures are located in the outer layer of 



Particle Acceleration in Reconnection Outflows 



11 



(a) AIA 131 A. 04:44:47 



(b) AIA 131 A. 04:48:59 (c) RHESSI 04:46:50-04:51:10 



X-roy Contours 04:42:30-04:46:50 04:46:50-04:51:10 
6-8 keV, Green 33, 95% 3-6 keV 

1 2 - 1 5 keV, Blue Doited 43, 95% 12-15 keV 





106 

104 

102 

100 

98 

96 

94 
59 
58 
57 
56 
55 



A 04:46:50-04:51:10 
+ 04:42:30-04:46:50 
o 04:38:10-04:42:30 
X 04:33:50-04:38:10 



220 240 260 
-Solor Y (orcsec) 



220 240 260 
-Solor Y (orcsec) 



220 240 260 280 
-Solor Y (orcsec) 




E (keV) 



Fig. 10. — Energy dependence of the double coronal X-ray sources likely located in the oppositely directed reconnection outflows, (a) 
and (b) RHESSI contours overlaid on AIA 131 A images. The higher energy sources (blue dotted contours at 43% and 95% of the image 
maximum) are closer toward one another than the lower energy sources (green contours), (c) Contoured 6-8 keV image at the time of (b). 
(d) and (e) Projected heights along the flducial Cut of the centroids of the two sources enclosed by the white dashed lines in (c). The 
result is color-coded for four consecutive times overlaid with linear flts in dashed line. The red horizontal error bars indicate energy bins. 



the flare arcade and the lower portion of the contracting 
cusps (ray-like features) above it. The EM in the latter 
is nearly four orders of magnitude lower than in the cen- 
tral arcade, similar to those of hot fans of rays observed 
by Yohkoh (|Svestka et alJ |1998[ ) . The considerable in- 
crease of EM and thus density in the flare loops from 
panel (d) to (e), prior to the impulsive phase, suggests 
chro mospheric evaporation driven by thermal conduc- 
tion (iLiu et alJl2QQ9bl:lBattaglia et alJ[2QQ9[ ) or by Alfven 
waves (|Haerendelll2009[ r rather than electron beam heat- 
ing that may be important later in the presence of HXR 
footpoints. 

To better follow the history of T and EM, we obtained 
their space-time plots, as shown in Figures [TTlf g) and 
(k), from corresponding maps using Cut 0. We iden- 
tified the maximum temperature and EM at each time 
and show their projected heights as orange and black 
symbols in panel (p). In general, both peaks follow the 
same upward- downward- upward motions as the loop-top 
emission centroids. The temperature peaks are near the 
RHESSI 16-25 keV centroids (blue), while the EM peaks 
are close to or slightly lower than the AIA 131 A centroids 
(purple) at lower heights. 

The peak offset of the temperature and EM can be 
better seen in their height distributions at selected times 
shown in Figures [TTlf h) and (1). It explains the observed 
energy dispersion of the X-ray loop-top centroids shown 
in Figure [11(e) because of the exponential shape of ther- 
mal bremsstrahlung spectra (equation (|4])). As shown 
by the green and red lines in the middle of Figure [9l 
a higher temperature but lower EM produces a harder 
(shallower) spectrum of lower normalization that dom- 
inates at high energies, while a lower temperature but 
higher EM produces a softer (steeper) spectrum of higher 
normalization that dominates at low energies. The tem- 
perature peak being located above the EM peak thus 
shifts the higher energy X-rays toward greater heights. 
Otherwise, if the temperature and EM peaks are cospa- 
tial, they would dominate X-rays at all energies and there 
would be no separation of centroids with energy. We note 
in Figures \TT\ h.) that the EM decreases more gradually 



near the flux rope. This may lead to the less pronounced 
energy dispersion of the upper coronal source there than 
the lower (loop-top) source (see Figure [TQ|) . 

As a proof of concept, we modeled the observed X- 
rays of energy E with thermal bremsstrahlung radiation 
(jTandb erg-Hanssen & Emslie 1988, p. 114): 



/sxR (X (EM) 



exp{-E/kT) 

eVt 



g{E/kT), 



(4) 



where EM is the emission measure and g{E/kT) = 
{kT/Ef^^ is the Gaunt factor. The resulting X -ray pro- 
files at two selected times are shown in Figures fTTTi) and 
(m) for the corresponding T and EM profiles. As ex- 
pected, the higher energy 20 keV emission (blue) is domi- 
nated by the temperature peak at a higher altitude, while 
the 3 keV emission (red) is dominated by the EM peak at 
a lower altitude. Their emission peaks, marked by open 
circles, differ in height by A/i = 1" and AQ" for the two 
times. Their peak heights are repeated in panel (p) and 
are close to those of observed loop-top centroids. 

As shown in Figures fTTTg) and (p), the high temper- 
ature region and particularly the temperature peak are 
close to the X-ray loop-top centroids but always below 
the reconnection site (green dashed line). This indicates 
that primary plasma heating takes place in reconnection 
outflows. To illustrate this, we can follow a contracting 
loop and the temperature variation it senses as it trav- 
els away from the reconnection site. A selected track 
from Figure [71(c) is shown here as the blue dotted arrow. 
The temperature history on its path, as shown in Fig- 
ures [Tlfh) and (j), indicates rapid heating at an average 
rate of 3.2MKmin~^ from 3.5 to 28 MK over 8 minutes 
and A/i = 90 Mm. This example gives us a sense of the 
heating rate averaged along the line of sight, not neces- 
sarily of a specific loop. 

Likewise, the track of a slow loop shrinkage at 
lower altitudes reveals cooling at an average rate of 
-0.32MKmin"^ from 25 to 7.9 MK within 54 minutes. 
The rate is — IMKmin"^ earlier during the impulsive 
phase in the 28-9 MK range and — O.lMKmin"^ later 
during 07:00-08:30 UT in the 14-7 MK range. These 
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Fig. 11. — Maps of temperature T (a-c) and emission measure EM (d-f) inferred from AIA filter ratios at three selected times on 
logarithmic scales, (g) and (k) Space-time plots of T and EM maps from Cut as shown in (b). The white dotted lines are temperature 
and emission measure from GOES data as shown in Figures [Tj^d) and (e). The first side panels ((h) and (1)) on the right show the height 
distributions of T (red, top scale) and EM (black, bottom scale) at the time marked by the triangle on the left. The second side panels 
((i) and (m)) show the modeled corresponding thermal X-ray emission at two photon energies, 3 and 20 keV. The open circles mark the 
heights of the emission peaks that are also shown in (p). (j) and (o) Temporal profiles T and EM, obtained from a horizontal slice of the 
corresponding space-time plots above at height /iref = 60^^ marked by the cross. The arrowed blue and light blue dotted lines in (g) are 
selected tracks of a fast contraction and slow shrinkage from Figures [Jl and [6] The temperatures sampled by these tracks are shown as a 
function of time and distance in (j) and (h), respectively, indicating heating and cooling experienced by these loops as they travel, (p) 
Projected heights of the maxima of T (orange) and EM (black) at each time in the corresponding space-time plots. Overlaid here and in 
the above space-time plots are the heights of the 16-25 keV (blue) and 131 A (purple) loop-top centroids in small symbols from Figure[4{e), 
and of the inferred reconnection site in green dashed line from Figure [Tfc). The vertical dashed lines mark the impulsive phase. 

2. Originating from the inferred magnetic reconnec- 
tion site within the current sheet are bi-directional 
outflows in the forms of plasmoids and cusp- 
shaped loops. The upward ejections have a me- 
dian initial velocity of 320kms~^ and a maximum 
of 1050 kms~^, while the concurrent downward 
cusp contractions are ~2/3 slower with a median 
of — 210kms~^. Even faster contractions up to 
— 920kms~^ occur after the flux rope eruption for 
another 10 h, and a sizable fraction of 21% of them 
have speeds >4 00kms~^, twice faster tha n previ- 
ously reported ([Savage fc McKenzid l201ll ). Such 
high velocities are comparable to expected coronal 
Alfven speeds of ~1000kms~^ and sound speeds 
of 520-910 kms-^ for a 10-30 MK flaring plasma. 
At lower altitudes, flare loops persistently shrink 



cooling rates are somewhat lower than previous ly found 
at even lower temperatures (jVrsnak et al.ll2006l ). 

6. CONCLUSION 

6.1. Summary 

We have presented detailed EUV and X-ray observa- 
tions of the 2012 July 19 M7.7 flare, focusing on the sig- 
natures of magnetic reconnection and associated energy 
release. We summarize our findings and discuss their 
implications as follows. 

1. The V-shaped EUV emission on the trailing edge 
of the flux rope CME and the underlying, inverted 
V-shaped flare loops, both associated with distinct 
X-ray emission, suggest two oppositely oriented Y- 
type null points with a vertical current sheet formed 
in between (Figures [2] and [5j) . 
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at typical initial velocity of — 17kms~^ and gradu- 
ally decelerate to about —5 kms~^ within 0.5-2.5 h. 
They are all evidence of reconnection outflows dur- 
ing different stages (Figures [6] and [7]). 

3. The double coronal X-ray sources are spatially sep- 
arated and located in the regions of bi-directional 
plasma outflows. The highest temperature is found 
near the loop-top X-ray source well below the re- 
connection site. This suggests that primary plasma 
heating and particle acceleration take p l ace i n 
the reconnection outflow regions dHolmadEoil), 
rather than at the reconnection site itself. Models 
with this ingredient w ere proposed long ago (e.g., 
iForbes fc Priestlll983f ). but solid observational ev- 
idence as presented here has been lacking (Fig- 
ures [TO] and [11]) . 

4. An energy dispersion is present in the loop-top po- 
sition with higher-energy X-ray and hotter EUV 
emission located at greater heights over a large 
range of Ah > 26 Mm. The 25-50 keV nonther- 
mal emission lies 15 Mm (twice that of the Masuda 
flare) above the 6-10 keV thermal emission. This 
agrees with the expected trend of softer electron 
spectra in the nonthermal regime and lower tem- 
peratures in the thermal regime being associated 
with earlier energized loops, which are located fur- 
ther below the primary locus of energy release. The 
upper coronal X-ray source has an opposite trend 
because it is located in the oppositely directed re- 
connection outflow (Figures [3l [U and [TO]). 

5. Prior to the recently recognized descent followed by 
a continuous ascent, the overall loop-top emission 
experiences an initial ascent for nearly an hour. 
This is the first time that such motions, including 
the descent, are observed simultaneously from EUV 
to HXRs covering a wide range of temperatures of 
1-30 MK. The transition from ascent to descent co- 
incides with the rapid acceleration of the flux rope 
CME. The loop-top descends at ~10kms~^, about 
50% slower in EUV than in X-ray (Figure H]). 

6. The flare impulsive phase starts when a rapid ve- 
locity increases occur for the overall loop-top de- 
scent, the individual loop shrinkages, and the up- 
ward plasmoid ejections. This is delayed by 10 
minutes from the initial loop-top descent, imply- 
ing that the energy release rate is more intimately 
correlated with these velocities than the loop-top 
position (Figure [4]). 

6.2. Proposed Physical Picture 

We propose the following physical picture to tie to- 
gether the observations. The major points, including 
the location of particle acceleration and the three-stage 
(up-down-up) motion of the loop-top X-ray source, are 
sketched in Figure [12] 

The earlier confined C4.5 flare leads to the formation 
of a flux rope and cusp-shaped loops underneath it. Af- 
ter six hours of sl ow evolution, the flux ro pe becomes 
unstable and rises (jPatsourakos et al]l2013l ). A vertical 
current sheet forms between two Y-type null points at 
the lower tip of the flux rope and the upper tip of the 



underlying cusp. Magnetic reconnection ensues within 
the current sheet, leading to the eruptive M7.7 flare. 

Magnetic reconnection produces bi-directional outflows 
in forms of the observed plasmoids and contracting loops. 
Plasmoids ar e flux tubes formed by the tearing mode 
(|Furth et al.l [T963t) with a guiding fleld along the cur- 
rent sheet. They are magnetic islands in two dimension 
or when the current sheet is seen edge-on. These out- 
flows are driven by the magnetic tension force of the 
highly bent, newly reconnected fleld lines, as in those 
pointed cusps. The outflows generally decelerate, as ob- 
served here, when they run into the ambient corona and 
when the contracting loops relax to less bent shapes with 
reduced magnetic tension. This implies that the low- 
altitude slow shrinkages could be the late stages of de- 
celerated high-altitude fast contractions. 

Several mechanisms can operate in the reconnection 
outflows and contribute to particle acceleration and 
plasma heating. Turbulence or plasma waves^ for ex- 
ample, can be generated by the interaction of the high 
speed flows with the ambient corona. Upon cascading 
to smaller scales (comparable to the gyro-radii of back- 
ground particles) at some distance from the reconnection 
site, tur bulence can accelerate the pa r ticles and heat the 
plasma ("Hamilton & Petrosian"1992'; 'Miller et al."1996l: 
Chandran 2003; Petrosian & Liu 2004; Petrosian e t al.l 
I2006i : .Jiang et al.i .2009» : .Fleishman Toptvgin .20131 ). 
Additional contribution to particle acceleration and/or 
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gas dynamic s hocks within contracting flux tubes 
(jLongcope et al.l [2009). and the first- order Fermi 
and betatron mechanisms with in the collapsing traps 
formed by contracting l oops ([Somov fc Kosugi 1997; 
Karhckv &: IWdl [2003 : iKarlickvl 120061 : IGradv et al.l 



As the flux rope evolves into its fast rise and eruption 
stage, the current sheet could become suffl ciently long 
causing signiflcant tearing instability (e.g., iBarta et"al] 
I2OO8I ). and/or become thinner than the ion skin depth 
leading to a transition from Sw eet-Parker reconnec tion to 
collisionless Hall reconnection (jCassak et al.|[2QQ6[ ) . Both 
can result in enhanced rates of reconnection and energy 
release. The outflow velocity will also increase, produc- 
ing stronger heating and particle acceleration by one or a 
combination of the above three mechanisms (turbulence, 
shock, and betatron) whose efficiency is expected to be 
positively correlated with the outflow velocity. For exam- 
ple, faster outflows can generate stronger turbu lence that 
can l ead to stronger particle acceleration (see iPetrosianI 
[20T2I . for a comparison between acceleration by turbu- 
lence and shocks). This explains the observed temporal 
correlation of the increase in outflow velocity and the 
onset of the impulse phase and HXR burst. 

The loop-top source seen in X-ray and EUV represents 
the collective emission of the ensemble of flare loops, 
each of which undergoes post-reconnection contraction or 
shrinkage. The loop- top emission centroid is the average 
position of this ensemble, which in the thermal regime 
is determined by the convolution of the spatial distribu- 
tions of temperature and emission measure. Such dis- 
tributions are influenced by several competing processes, 
including energization of new loops produced by the up- 
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Fig. 12. — (a) Schematic of the proposed flare model in which particle acceleration and plasma heating take place in the reconnect ion 
outflow regions away from the reconnection site, (b)-(e) Temporal development of reconnection and the up-down-up three-stage motion of 
the loop-top (LT) X-ray source. We show only the portions of reconnected magnetic fleld lines near the reconnection site, above and below 
which each pair of fleld lines are identifled by the same color during their relaxations. Note that the higher contraction speeds (e.g., of the 
green loop) from (c) to (d) is associated with the loop-top source descent. The upper coronal source is not hatched during the impulsive 
and decay phases, indicating its non-detection due to weak emission. We assume that bi-directional reconnection outflows and relaxations 
of fleld lines persist throughout the event, but the upward component is beyond the AIA FOV during the late phase and thus not observed. 



ward developing reconnection, and downward contrac- 
tions and cooling of previously reconnected loops. We 
suggest that the interplay of these processes results in the 
up-down-up migration of the loop-top source. Around 
the early impulsive phase, the observed higher velocities 
of loop contractions tend to shift the overall loop-top 
source downward, as depicted in Figures [T2fc) and (d). 
Within the contracting loops, conductive cooling can be 
considerably reduced by stron ger turbulence produced by 
faster reconnectio n outflows (jChandran fc CowlevUlQQSl : 
iJiang et aI]l2QQ6f ). Slower cooling helps the X-ray /EUV 
emission of these hot loops last longer during their con- 
tractions and further contribute to the loop-top descent. 
Before and after the early impulsive phase, the situation 
could be different. Loop contractions are slower and asso- 
ciated faster cooling can make these loops cool below the 
instrumental temperature passband more rapidly. The 
upward development of reconnection could thus domi- 
nate, leading to the upward loop-top migration. This in- 
terpretation is supported by the observed temporal coin- 
cidence of the highest velocities of individual loop shrink- 
ages and of the loop-top descent at the impulsive phase 
onset. Like cooling, chromospheric evaporation can also 
play a role. After the initial chromospheric evaporation 
associated with energization of newly reconnected loops, 
if it continues operating during the subsequent loop con- 
tractions, the densities and thus emission measures of 
these loops would keep rising and contribute to the down- 
ward loop- top drift. Otherwise, if it rapidly subsides, 
it would assist the upward loop-top drift. Investigating 
such a complex interplay would require detailed numeri- 
cal modeling, which is beyond our scope here. 

6.3. Discussion 

The fast downward contractions of EUV loops are 
best seen during the decay phase, likely because of their 
greater distances traveled making them easier to be de- 
tected. Likewise, different environments may explain the 
factor of two difference in the median deceleration be- 



tween the early and later fast contractions (Figures [71^ f) 
and (g)). The earlier, compact flare loops of stronger 
magnetic field can produce stronger resista nce (e.g., by 
the L orentz force from a reverse current; iBarta et al.l 
l2QQ8f ) to quickly brake the contracting loops impinging 
on them from above, while the later, large flare arcade 
of weaker field at greater heights and greater travel dis- 
tances may allow more gradual deceleration. 

The X-ray spectra are essentially thermal during the 
decay phase. This suggests that the late phase contrac- 
tions are associated with plasma heating rather than par- 
ticle acceleration. This is expected because the mag- 
netic field strength decreases with height, and so does 
the available magnetic energy to be released by reconnec- 
tion. Such prolonged heating may contribute to the well- 
known slowe r than expected cooling of flare plasma ( see 
Figure [U^d) ; iMcTiernan et al.l 119931 : iJiang et al.l l2QQ6f ) . 

The fast high-altitude contractions and slow low- 
altitude shrinkages appear as two statistically distinct 
populations. The median final velocity of the former is 
more than four times the median initial velocity of the 
latter (—73 vs. — 17kms~^; Figures [6] and [7|) • We have 
interpreted the latter as the late stages of the former 
that have decelerated, but we rarely see a continuous 
track decelerating from >400kms~^ to a few kms~^. 
Alternatively, this distinction may be due to inho mo- 
geneity of reconnection (see, e.g., lAsai et aTl 1200 4^ in 
the current sheet above the flare arcade that is primar- 
ily oriented along the line of sight. The persistent slow 
shrinkages could result from reconnection at a moder- 
ate rate throughout the current sheet, which produces 
flare loops forming the arcade and double ribbons at 
their footpoints. The episodic fast contractions could 
be signatures of plasmoids within the current sheet, as- 
sociated with enhanced reconnection and energy release 
rates. These fast contracting loops are dispersed along 
the line of sight amidst slowly shrinking loops, producing 
the observed effects (e.g.. Figure (Tfb)). 

There are other physically different (though not nec- 
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essarily unrelated) phenomena that share similar obser- 
vational signatures and should not be confused with the 
contractions of newly reconnected loops studied here. For 
example, pre-existing coronal loops can co ntract at typ- 
ical speeds < lQQkms~ ^ during eruptio ns ()Liu fc Wand 
I2Q1QI : iLiu et aL 2Q12d: iSun etaIll2Q12p. int erpreted as 
implosion ()Hudsonll2QQQl : iJanse fc Lowll2QQ7[ ) due to the 
rapid release of magnetic energy and thus reduction of 
magnetic pressure within the eruption volume. When 
a CME eruption initiates at an elevated height, it can 
produces a do wnward push to displace ambient coro- 
nal loops (e.g.. iLiu et al.|[2Q12d see their Figure 12(h), 
at — 60kms~^). Even slower shrinkage at — 3kms~^ 
can occur in active region loops due to gradual cooling 
(|Wang et al.lll997l ). Sometimes, conjugate X-ray foot- 
points approach each other while the loop-top source de - 
scends (jJi et al.ll2QQ6l : lLiu et al.ll2QQ9aHYang et al.ll2QQ9[ ). 
likely because of reconnection progressing toward less 
sheared loops at lower altitudes. We have not found 
these signatures in the flare under study here, except 
for a few overlying existing loops that contract at about 
— 150kms~^ during 05:40-05:50 UT near the end of the 
impulsive phase. In addition, some high-lying loops 
continue to expand and erupt at ~100kms~^ up to 
08:00 UT, three hours after the original flux rope erup- 
tion. Examination of such features and quantitative com- 
parison of these observations with particle acceleration 
models will be subjects of future investigations. 
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Hudson, Sasha Kosovichev, Judy Karpen, Rick DeVore, 
Greg Slater, Tom Berger, and Mark Cheung for useful 
discussions. 

APPENDIX A 
OBSERVATIONS OF ALL AIA CHANNELS 

In the main text, we focused on the 131 A channel. For 
a complete temperature coverage, we examine all AIA 
channels here. As shown in Figure [T3l (top), each chan- 
nel has a generally broad response with one to a few 
peaks. In the order of approximately decreasing temper- 
ature response are 193, 131, 94, 304, 335, 211, and 171 A 
channels. 

Figure [13] shows space- time plots of all EUV channels. 
On the left are base ratio from every other frame of short 
exposure so that bright flare loops are not over-exposed, 
while on the right are running ratio from frames of regu- 
lar exposure that is needed for detecting faint emission, 
including fast loop contractions, above the flare arcade. 
On the left, the 131 and 94 A channels capture the hot, 
diffuse loop-top emission, while the flare arcade under- 
neath it is best seen in cooler channels. Both regions ap- 
pear in the 193 A channel because of its dual response to 
hot Fe XXIV emissi on peaking at 18 MK and cool Fe XII 
emission at 1.6 MK (jO'Dwver et al.ll2010[ ). On the right- 
hand side, downward fast contractions and slow shrink- 



ages can be identified in the three hottest channels. In 
cooler channels, cooling condensations rain down the ar- 
cade loops since 06:00 UT at up to ~170kms~^ or 90% 
of the free fall speed, faster than typi cal coronal rain 
(jAntolin fc Rouppe van der Voortl[2012l ). 

In Figure [T3l (left), the erupting flux rope is best seen 
in the hot 131 and 94 A channels as a bright, accelerat- 
ing track, indicating its high temperature, as also shown 
in Figure [TTJ Its evolution from gradual rise to impul- 
sive eruption is accompanied by the rapid increase in 
X-ray flux and the loop-top's ascent-to-descent transi- 
tion. This is in line with obser ved synchronous eruption - 
accelerations and flare onsets (IZhang et aI]l200lL 120021 : 
iTemmer et alll2008l : iLiii et al.ll2(nT^ 

Figures [U^b) andlH^a) show temporal emission profiles 
taken from the base-ratio space-time plots at /iref = 60'' 
marked by a cross. Since 04:17 UT, the 131 and 94 A 
emission increases considerably, while other emission re- 
mains essentially flat. This indicates gentle heating early 
in the flare, as seen in the temperature history at this po- 
sition shown in Figure fTTT j). The dip in all channels at 
the onset of the impulsive phase (vertical dashed line) 
is due to the descent of the loop-top source. The large 
hump afterwards results from the growth of the flare ar- 
cade toward greater heights. This hump generally pro- 
gresses toward cooler channels, indicating a cooling se- 
quence, similar to those recently observed in flares, AR 
loops, and prominence condensations ( Woods et al. 2QTTI: 
IViall fc Klimchug [20T1 ILiu et al.l l2012bl : iBerger et al.l 
l2012f ). The double humps at 193 A are again due to 
its dual temperature response. 

Like 193 A, the 304 A channel responds to both hot 
Ca XVIII emission at 7.1 MK and cool He II emission 
at 0.05 MK. The dual response of these two channels 
leads to their unusually large descent of about 40% of 
the loop-top centroids (see Table [2]). This is because 
their maximum loop-top heights are dominated by hot 
emission, while the minimum heights are dominated by 
cool emission at even lower heights (Figure |4]). 

APPENDIX B 
DISAMBIGUATION WITH STEREO OBSERVATIONS 

The limb view of SDO/AIA leaves an alternative pos- 
sibility of the loop-top descent as the progression of en- 
ergization from high- to low-lying loops along the arcade 
in the east-west direction. This possibility is not sup- 
ported by STEREO-A (STA) that was 121° ahead of 
the Earth and observed the flare on the disk. When 
the X-ray loop-top descent occurs, the STA 195 A im- 
age (Figure [MT d)) shows a compact cusp-shaped loop, 
while the arcade system develops later (panels (e) and 
(f)). One could further argue that a temperature effect, 
i.e., cooling or heating, can make a progression appear 
as a compact source rather than an arcade. If so, the 
loops being cooled or heated would have appeared at the 
same location and been captured by AIA in cooler or 
hotter EUV channels or by RHESSI in X-rays. Instead, 
all available EUV and X-ray data covering a wide tem- 
perature range of 10^-10^ K show a consistent loop-top 
descent, which must be a real mass motion. 
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Fig. 14. Comparison of SDO/AIA 193 A (top) and STEREO-A 195 A (bottom) images of the flare. The STEREO images are scaled 
as if they are viewed from the same hehocentric distance as SDO. The hehographic grids are spaced at 5°. Note the compact brightening 
at the impulsive phase onset (pointed by the arrow), in contrast to the extended arcade system later. 
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